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Abstract: 

We have calculated the complete electroweak 0{a) radiative corrections to the Higgs- 
boson production process e"''e~ — > ttH in the electroweak Standard Model. Initial-state 
radiation beyond 0{a) is included in the structure-function approach. The calculation of 
the corrections is briefly described, and numerical results are presented for the total cross 
section. Both the photonic and the genuine weak corrections reach the order of about 
10% or even more and show a non-trivial dependence on the Higgs-boson mass and on the 
scattering energy. We compare our results with two previous calculations that obtained 
differing results at high energies. 
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1 Introduction 



In the electroweak Standard Model (SM) all fermions / receive their masses ruf via 
Yukawa couplings to the Higgs field. Splitting the Higgs field into its vacuum-expectation 
value V = (v^G^)^^/^ ~ 246 GeV and the physical excitation H(x), the Yukawa term in 
the SM Lagrangian reads Lyuk = — Hf"mfi)- + H{x)/v)iljfipf. Thus, the Yukawa coupling 
strength is predicted to be m^/f at tree level, and the experimental determination of the 
Higgs-fermion couplings represents an important test of the mass generation via the Higgs 
mechanism. Since the top quark is the heaviest of all fermions, it is supposed to play a key 
role in a theory of fermion masses. Therefore, the measurement of the top-quark Yukawa 
coupling is of particular interest. For not too large Higgs-boson masses, Mh ~ 100- 
200 GeV, a promising process for this task is e+e" — > ttH, as already pointed out in 
Ref. [p. To achieve a measurement with a precision of the order of ~ 5%, however, an 
e"'"e~ linear collider (LC) with high centre-of-mass (CM) energies (i/i ~ 800-1000 GeV) 
and high luminosity (L ~ 1000 fb^^) is required [|2]. A better accuracy might be obtained 
by a simultaneous fit of various Higgs-boson parameters to the whole profile of the Higgs 
boson at a LC [E]. Moreover, an investigation of the process e"'"e~ — > ttH might be useful 
for setting bounds on non-standard physics [|!T||!lj in the top-quark-Higgs-boson coupling. 

A determination of the top-quark Yukawa coupling at the level of a few per cent 
requires both a proper understanding of the background [ |S] to the decaying ttH final 
state and a theoretical prediction of the e"'"e~ ttH signal cross section within per-cent 
accuracy. Thus, radiative corrections have to be controlled within this accuracy, a task 
that is rather complicated for a process with three massive unstable particles in the final 
state. As a first step, the process e"^e~ — >■ ttH can be treated in the approximation of 
stable top quarks and Higgs bosons in the final state. The corresponding lowest-order 
predictions are already known for a long time [E]. These results were supplemented by 
the 0{as) QCD corrections to the total production cross section in the SM, first within the 
"effective Higgs-boson approximation" [ 7^ that is valid only for small Higgs-boson masses 
and very high energies, and subsequently [^^^ based on the full set of QCD diagrams in 
0{as)- The QCD corrections to Higgs-boson production in association with heavy-quark 
pairs (tt/bb) in the minimal supersymmetric SM were discussed in Refs. [ ITmillllT^ .^ 
Besides corrections to total cross sections, the QCD corrections to the Higgs-boson energy 
distribution were discussed in Ref. [E!, both for the SM and its minimal supersymmetric 
extension. Recently first results for the electroweak 0{a) corrections to e"'"e~ ttH in 
the SM have been presented in Refs. [|T3j and [CH]. These calculations were found to 
agree for small energies but to differ at high energies. 

In this paper we present results of a further, completely independent calculation of the 
0{a) electroweak corrections, which is additionally improved by the leading higher-order 
corrections from initial-state radiation. Details on this calculation, which we have also 
implemented in a Monte Carlo event generator, will be given elsewhere. Here we sketch 

^In Ref. [E] only the QCD corrections to the photon-exchange channel for ttH production were taken 
into account. 

^In Ref. only the photon-exchange channel is corrected, so that the relative correction given there 
coincides with the one obtained in Ref. [ 9 where the SM process is treated analogously. The calculation 
of Ref. includes the full set of 0{as) QCD diagrams, and in Ref. [ 12 the SUSY-QCD corrections 
to e+e~ ttH are considered. 
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only the main ingredients. Moreover, we compare our results with those of Refs. [ El and 
[ I15j . While we find good agreement with the results of the recent calculation of Ref. [ IT^ 
at all considered energies, our results differ from those of Ref. [El at high energies and 
close to threshold. Moreover, we could reproduce the QCD corrections of Ref. [[S] from 
our results on photonic final-state radiation within statistical integration errors. 

2 Method of calculation 

We have calculated the complete 0{a) electroweak virtual and real photonic correc- 
tions to the process e"'"e~ —>■ ttH, following the same strategy as used in our previous 
calculation of the corrections to e"'"e~ uuH 

The calculation of the one-loop diagrams has been performed in the 't Hooft-Feynman 
gauge both in the conventional and in the background-field formalism using the conven- 
tions of Refs. [Ej and [El, respectively. The renormalization is carried out in the on-shell 
renormalization scheme, as described there. The electron mass mg is neglected whenever 
possible. 

The calculation of the Feynman diagrams has been performed in two completely inde- 
pendent ways, leading to two independent computer codes for the numerical evaluation. 
Both calculations are based on the methods described in Ref. [ Ej ■ Apart from the 5-point 
functions the tensor coefficients of the one-loop integrals are recursively reduced to scalar 
integrals with the Passarino-Veltman algorithm [E] at the numerical level. The scalar 
integrals are evaluated using the methods and results of Refs. [ [d ED] , where ultraviolet 
divergences are regulated dimensionally and IR divergences with an infinitesimal photon 
mass. The 5-point functions are reduced to 4-point functions following Ref. [1211, where a 
method for a direct reduction is described that avoids leading inverse Gram determinants 
which potentially cause numerical instabilities. The two calculations differ in the following 
points. In the first calculation, the Feynman graphs are generated with FeynArts version 
1.0 [[23- Using Mathematica the amplitudes are expressed in terms of standard matrix 
elements and coefficients of tensor integrals. The whole calculation has been carried out 
in the conventional and in the background-field formalism. The second calculation has 
been done with the help of FeynArts version 3 [ , and the analytical expressions have 
been generated by FormCalc [|21I and translated into a C code. The scalar and tensor 
coefficients, in particular those for the 5-point functions, have been evaluated by own 
routines, as described above. 

The results of the two different codes, and also those obtained within the conventional 
and background-field formalism, are in good numerical agreement (typically within at least 
12 digits for non-exceptional phase-space points and double precision) . The agreement of 
the results in the conventional and background-field formalism, in particular, checks the 
gauge independence of our results. 

The matrix elements for the real photonic corrections are evaluated using the Weyl- 
van der Waerden spinor technique as formulated in Ref. [ '25^ and have been successfully 
checked against the result obtained with the package Madgraph [ 26 . The soft and 
collinear singularities are treated in the dipole subtraction method [l27j. Beyond 0{a) 
initial-state-radiation (ISR) corrections are included at the leading-logarithmic level using 
the structure functions given in Ref. [|2H] (for the original papers see references therein). 
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The phase-space integration is performed with Monte Carlo techniques in both com- 
puter codes. The first code employs a multi-channel Monte Carlo generator similar to 
the one implemented in RacoonWW [123 EO] and Lusifer [EI], the second one uses the 
adaptive multi-dimensional integration program VEGAS [E2]- 

3 Numerical results 

For the numerical evaluation we use the following set of SM parameters [ j33j , 





= 1.16639 X 10-^GeV"^ 


a(0) 


= 1/137.03599976, 


as(Mz) 


= 0.1172, 


Mw 


= 80.423 GeV, 


Mz 


= 91.1876 GeV, 






me 


= 0.510998902 MeV, 




= 105.658357 MeV, 


m^ 


= 1.77699 GeV, 


mu 


= 66 MeV, 


rric 


= 1.2 GeV, 


mt 


= 174.3 GeV, 


md 


= 66 MeV, 


ms 


= 150 MeV, 


mb 


= 4.3 GeV. 



(3.1) 

We do not calculate the W-boson mass from but use its experimental value as input. 
The masses of the light quarks are adjusted to reproduce the hadronic contribution to the 
photonic vacuum polarization of Ref. [[34]. We parametrize the lowest-order cross section 
with the Fermi constant G^ (G^-scheme), i.e. we derive the electromagnetic coupling 
a according to ac^ = V^G^M^sl^/n. This, in particular, absorbs the running of the 
electromagnetic coupling a(Q^) from Q = to the electroweak scale {Q ~ Mz) into the 
lowest-order cross section so that the results are practically independent of the masses of 
the light quarks. In the relative radiative corrections, we use, however, a(0) as couphng 
parameter, which is the correct effective coupling for real photon emission. 

In the following we separate the photonic corrections, which comprise loop diagrams 
with virtual photon exchange in the loop and the corresponding parts of the counter terms 
as well as real photon emission, from the full electroweak 0{a) corrections; the remaining 
non-photonic electroweak corrections are called weak. Since the lowest-order diagrams 
involve only neutral-current couplings, but no W-boson exchange, this splitting is gauge 
invariant. Moreover, we separately discuss the higher-order [i.e. beyond 0(a)] ISR cor- 
rections that are obtained from the convolution of the lowest-order cross section with the 
leading- logarithmic structure functions, but with the 0{a) contribution subtracted. The 
higher-order ISR is included in the electroweak corrections shown in the following plots. 

In addition to the electroweak corrections, we also include results on the QCD correc- 
tions which can be deduced from the part of the photonic corrections that is proportional 
to Qt, where Qt = 2/3 is the relative electric charge of the top quark, i.e. from the 
final-state radiation (FSR). The QCD correction is obtained from these corrections upon 
replacing the factor Q^a by CpOs = 4as(/i^)/3. Following Ref. [|Hj, the QCD renormaliza- 
tion scale /z is set to the CM energy in the following, and the running of the strong coupling 
is evaluated at the two-loop level (MS scheme) with five active flavours, normalized by 
as(M|) as given in Eq. ()3.1|) . For y/s = 500,800, and 1000 GeV the resulting values for 
the strong coupling are given by as(M|) = 0.09349,0.08857, and 0.08642, respectively. 

In this letter, we consider merely total cross sections without any cuts; distributions 
will be discussed elsewhere. For reference we give some numbers for the total cross section 
in lowest order, o"tree, the cross section including electroweak and QCD corrections, o"corr, 
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Figure 1: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-sclieme for a CM energy ^/s = 500 GeV 



and the various contributions to the relative corrections defined as 5 = a/atiee — 1 in 
Tabled The last numbers in parentheses correspond to the Monte Carlo integration 
error of the last given digits. In Figures we show the lowest-order cross section 
as well as the corresponding corrections as a function of the Higgs-boson mass for the 
typical LC CM energies y/s = 500 GeV, 800 GeV, and 1 TeV. In Figures HHi the results 
are illustrated for fixed Higgs-boson masses of Mh = 115 GeV, 150 GeV, and 200 GeV. 

Away from the kinematic threshold at y/s = 2mt + Mh the total cross section is 
typically of the order of 0.5-3 fb and becomes maximal in the energy range between 
700 GeV and 1 TeV for small Higgs-boson masses. As already discussed in Refs. [ IHj in 
detail, the QCD corrections are positive and rather large in the threshold region {^/s > 
2mt + Mh), where soft-gluon exchange between in the tt system leads to a Coulomb-like 
singularity. Away from threshold this singularity is diluted by the larger phase space, since 
the singularity demands a low relative velocity of the quarks. Averaging the singular factor 
over the phase space leads to the following threshold behaviour of the QCD correction [ 



5qcd - 5fsr ~ A = , (3.2) 

where A is the maximal quark velocity in the tt rest frame. The FSR corrections show 
the same behaviour^, however, suppressed by the factor Q^a/{4:as/3) ~ 0.02. Although 
the relative QCD correction becomes rather huge close to threshold, it should be realized 
that the total cross section decreases rapidly there. In the region above threshold, where 
the cross section is largest, the QCD correction is only of the order of a per cent. For 



^Because of the dominance of the ISR corrections this is not visible in the foUowing figures. 
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/TT n 
Vs [GeVJ 




500 




Mh [GeV] 


115 


125 


140 


O-tree [fb] 


0.47901(7) 


0.23150(3) 


0.038189(6) 


(Tcorr [fb] 


0.6506(6) 


0.3401(4) 


0.0713(1) 


5qcd[%] 


46.5(1) 


59.2(2) 


103.4(3) 


Vot[%] 


-30.64(1) 


-34.69(1) 


-44.51(1) 


'5hoISR[%] 


4.25(3) 


5.53(3) 


9.28(3) 


(^weak [%] 


15.70(3) 


16.85(3) 


18.51(3) 



x/i [GeV] 








800 






Mh [GeV] 


115 


150 


200 


250 


300 


350 


O-tree [&] 


2.7004(4) 


1.7406(3) 


0.9217(1) 


0.46076(6) 


0.20432(3) 


0.07165(1) 


(Tcorr [fb] 


2.541(1) 


1.6076(7) 


0.8443(4) 


0.4251(2) 


0.1904(1) 


0.06648(6) 


5qcd[%] 


-0.87(2) 


0.36(3) 


2.43(4) 


5.03(4) 


8.58(6) 


14.12(6) 


5phot[%] 


-5.30(1) 


-8.26(1) 


-12.54(1) 


-16.76(1) 


-21.54(1) 


-27.55(1) 


^hoISR[%] 


-0.41(3) 


-0.19(3) 


0.22(3) 


0.79(3) 


1.59(3) 


2.86(3) 


5weak[%] 


0.69(1) 


0.45(1) 


1.49(1) 


3.19(1) 


4.56(1) 


3.36(5) 



[GeV] 






1000 






Mh [GeV] 


115 


150 


200 


250 


300 


350 


CTtree [fb] 


2.2594(3) 


1.6208(2) 


1.0356(2) 


0.6643(1) 


0.41894(6) 


0.25510(4) 


CTcorr [fb] 


2.061(1) 


1.4311(7) 


0.8900(5) 


0.5639(3) 


0.3510(2) 


0.2045(2) 


5qcd[%] 


-4.76(2) 


-4.37(2) 


-3.72(3) 


-2.99(3) 


-2.26(4) 


-1.73(5) 


^phot[%] 


-0.34(1) 


-2.89(1) 


-6.35(1) 


-9.41(1) 


-12.43(1) 


-15.55(1) 


5hoISR[%] 


-0.61(3) 


-0.52(3) 


-0.34(3) 


-0.11(3) 


0.19(3) 


0.57(3) 


^weak[%] 


-3.06(2) 


-3.91(2) 


-3.66(3) 


-2.61(2) 


-1.71(2) 


-3.13(3) 



Table 1: Lowest-order cross section for e+e" — > ttH in the G^-scheme, Utree, cross section 
including electroweak and QCD corrections, o-com and various contributions to the relative 
corrections 6 (as described in the text) for various Higgs-boson masses at ^/s = 500 GeV, 
800 GeV, and 1 TeV 
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Figure 2: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-scheme for a CM energy y/s = 800 GeV 
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Figure 3: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-scheme for a CM energy ^/s = 1 TeV 
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Figure 4: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-scheme for a Higgs-boson mass Mh = 115 GeV 



a[fb] 



1 - 



0.1 - 



0.01 




tree 
corrected 



Mh = 150 GeV 



500 750 1000 1250 
G [GeV] 



(7tree 

40 



30 
20 
10 


-10 
-20 
-30 
-40 



1 ; 1 


1 1 

electroweak 




weak 




photonic - 


1 \ \ 


h.o. ISR 


1 \ V 

^ - - - ^ 


QCD ~ 






Mh = 150 GeV 


I' 1 


1 1 



1500 



500 750 1000 1250 
^ [GeV] 



1500 



Figure 5: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-scheme for a Higgs-boson mass Mh = 150 GeV 
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Figure 6: Lowest-order and corrected cross sections (l.h.s.) as well as relative corrections 
(r.h.s.) in the G^-scheme for a Higgs-boson mass Mh = 200 GeV 



CM energies and Higgs-boson masses far above threshold, the QCD corrections even turn 
negative and reduce the cross section by about 5% at ^/s = 1 TeV. This behaviour is 
expected from the effective Higgs-boson approximation [[7], as also argued in Ref. [fSl. 

The photonic corrections are negative for CM energies below 1 TeV and not too small 
Higgs-boson masses, i.e. not too far away from threshold. This is due to large virtual 
photonic corrections which are not cancelled by the corresponding real radiation owing 
to the decreasing phase-space volume for real-photon emission with increasing Mh or 
decreasing ^/s. Therefore, the photonic corrections are large and negative, in particular, 
at the relatively small scattering energy of y/s = 500 GeV, where they reach —25% to 
—65% for Mh = 100-150 GeV. In this situation resummation of the large photonic 
corrections is mandatory. The bulk of these contributions is included in our calculation 
(higher-order ISR) and for ^/s = 500 GeV amounts to 3-20%. Away from threshold, the 
photonic 0{a) corrections grow with increasing -y/i and reach +7% for Mh = 115 GeV 
and -y/s = 1.5 TeV. The higher-order ISR stays below 1% for y/s at least 150 GeV above 
threshold. 

The genuine weak corrections strongly depend on the scattering energy, while the 
dependence on the Higgs-boson mass is moderate. For y/s = 500 GeV they are about 15- 
20%, i.e. large and positive. For increasing CM energy they decrease and are at the per- 
cent level around ^/s ~ 800 GeV. For TeV energies they become more and more negative 
and reach the order of —10% around 1.5 TeV. While the weak corrections are smaller than 
the QCD corrections in the threshold region, both contributions are of similar size about 
250 GeV above threshold, i.e. near the peak of the cross section. At high energies the size of 
the weak corrections grows faster than the one of the QCD corrections. Such a behaviour 
is typical if Sudakov logarithms like aln^(s/M^) dominate the weak corrections at high 
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energies. Note that we use the G^-scheme, and that the weak corrections are shifted by 
~ ±10% when transformed to other schemes hke the a{0) or the a{s) schemes (see also 
Ref. [ilfij). The spikes at Mh = 2M-w,2Mz resuh from thresholds and are well-known 
from the process e"'"e~ ZH [I35j. 

Photonic and weak corrections partially cancel each other, and the resulting elec- 
troweak corrections increase very weakly with ^/s apart from the region close to thresh- 
old, where the electroweak corrections decrease fast with decreasing ^/s and reach of the 
order of —20% at threshold. Away from threshold, they are at the level of —5%, —8%, 
and —11% for Mh = 115 GeV, 150 GeV, and 200 GeV, respectively, and thus become 
increasingly negative with increasing Higgs-boson mass. 

4 Comparison with other calculations 

The results on the QCD corrections given in Tabled have been reproduced with the 
(publically available) computer code based on the calculation of Ref. [|Hj. We found 
agreement within the statistical integration errors. 

For a comparison of the electroweak 0{a) corrections with the results of Ref. [ IT^ 
we changed our input parameters to the ones quoted there and switched to the a(0)- 
scheme, where is ignored in the input and all couplings are deduced from a(0). In 
Table 121 we compare some representative numbers^ from the calculation of Ref. [^14j with 
the corresponding results from our Monte Carlo generator. The numbers in parentheses 
give the errors in the last digits of our calculation. The tree-level cross sections coincide 
within 0.03%. Most of the numbers for the one-loop corrected cross sections agree within 
1-2%, i.e. roughly within the estimated error of Ref. [El- However, for the corrected 
cross sections at ^/s = 2 TeV, i.e. at high energies, and the one very close to threshold, 
i.e. for i/i = 500 GeV and Mh = 150 GeV, we find differences of 4% and 7%, respectively. 
The same holds for the relative corrections. Ours are larger by about 4% at ^/s = 2 TeV 
and smaller by about 5% for the selected cross section close to threshold. 

Finally, we have also compared the electroweak 0{a) corrections with Ref. [El? where 
also the a(0)-scheme has been used. In Table E] we list the results of Table 2 of Ref. [ IT3] 
and the separate relative photonic and weak corrections of Table 3 of that paper together 
with the corresponding results from our Monte Carlo generator. Again the numbers in 
parentheses give the errors in the last digits. We reproduce the results for the lowest-order 
cross section within the integration errors, which are about 2-3 x 10~^. The results for 
the cross section including electroweak corrections coincide to better than 0.1% which is 
of the order of the integration error of the results of Ref. [ ■ The relative electroweak, 
weak, and photonic corrections agree also within 0.1%. This holds as well for the QCD 
corrections (not shown in Table E)). 

5 Summary 

We have presented results from a calculation of electroweak radiative corrections to 
the process e+e~ —>■ ttH, which is important for a precise determination of the top-Higgs 

^These numbers were kindly provided to us by Zhang Ren- You and You Yu quoting a statistical error 
below 1%. 
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Table 2: Total cross section in lowest order and including the full electroweak 0{a) 
corrections as well as the relative corrections for various CM energies and Higgs-boson 
masses for the input-parameter scheme of Ref. [E]- The statistical errors of Ref. [ Hi] 
are estimated by the authors to be below 1% (cf. footnote 4). 
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Mh 


O-troe [fb] 


O-O(a) [fb] 


5ew [%] 


^woak [%] 


5phot [%] 




[GeV] 


[GeV] 














600 


120 


1.7293(3) 


1.738(2) 


0.5 


16.5 


-16.0 


Ref. [Cg 




120 


1.7292(2) 


1.7368(6) 


0.44(3) 


16.49(3) 


-16.03(1) 


this work 


600 


180 


0.33714(4) 


0.3126(3) 


-7.3 


18.4 


-25.7 


Ref. [Cg 




180 


0.33714(5) 


0.3124(1) 


-7.34(3) 


18.38(3) 


-25.72(1) 


this work 


800 


120 


2.2724(5) 


2.362(4) 


3.9 


9.5 


-5.6 


Ref. [Cg 




120 


2.2723(3) 


2.3599(6) 


3.86(2) 


9.56(2) 


-5.70(1) 


this work 


800 


180 


1.0672(3) 


1.050(2) 


-1.6 


9.1 


-10.7 


Ref. [Cg 




180 


1.0668(2) 


1.0494(2) 


-1.63(2) 


9.18(2) 


-10.81(1) 


this work 


1000 


120 


1.9273(5) 


2.027(4) 


5.2 


5.8 


-0.6 


Ref. [cn] 




120 


1.9271(3) 


2.0252(5) 


5.09(2) 


5.78(2) 


-0.70(1) 


this work 


1000 


180 


1.1040(3) 


1.098(2) 


-0.5 


4.4 


-4.9 


Ref. [cn] 




180 


1.1039(2) 


1.0972(3) 


-0.61(2) 


4.41(2) 


-5.00(1) 


this work 



Table 3: Total cross section in lowest order and including the full electroweak 0{a) 
corrections as well as the relative electroweak, weak, and photonic corrections for various 
CM energies and Higgs-boson masses for the input-parameter scheme of Ref. [ [T3] 

Yukawa coupling. In detail, we have discussed the impact of photonic corrections at and 
beyond 0{a), the genuine weak 0{a) corrections, and the 0{as) QCD corrections. 

The photonic and weak corrections both reach the order of ~ 10% and show charac- 
teristic dependences on the Higgs-boson mass and on the scattering energy. Owing to a 
phase-space effect the (negative) photonic corrections reduce the cross section more and 
more when the threshold is approached, i.e. with increasing Higgs-boson masses. Close to 
threshold the resummation of the higher-order ISR corrections becomes mandatory. The 
weak corrections, which depend on the Higgs-boson masses only moderately, range from 
about +15% at a CM energy of 500 GeV to about —10% at 1.5 TeV. There are large 
cancellations between photonic and weak corrections, and the final size of the correc- 
tions depends strongly on the input-parameter scheme. These results clearly demonstrate 
the necessity to include the electroweak corrections in predictions adequate for a future 
high-luminosity e"'"e~ collider. 

We have compared our results with those of other groups. The QCD corrections have 
been successfully checked against previous calculations. The electroweak corrections agree 
well with the results of the recent calculation [E] but are at variance with the results of 
Ref. [mj at high energies and close to threshold. 
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